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Abstract 

We analyze the Hamilton- Jacobi action of gravity and matter in the limit where 
gravity is treated at the background field approximation. The motivation is to 
clarify when and how the solutions of the Wheeler-DeWitt equation lead to the 
Schrodinger equation in a given background. To this end, we determine when and 
how the total action, solution of the constraint equations of General Relativity, 
leads to the HJ action for matter in a given background. This is achieved by com- 
paring two neighboring solutions differing slightly in their matter energy content. 
To first order in the change of the 3-geometries, the change of the gravitational 
action equals the integral of the matter energy evaluated in the background ge- 
ometry. Higher order terms are governed by the "susceptibility" of the geometry. 
These classical properties also apply to quantum cosmology since the conditions 
which legitimize the use of WKB gravitational waves are concomitant with those 
governing the validity of the background field approximation. 
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1 Introduction 



There is a question that haunts all works in which both quantum matter effects 
and the dynamics of gravity play a crucial role, see e.g. P], E2L EjJ: What is the 
validity of the "semi-classical" Einstein's equations 

- g^R/2 = 8nG (*|T Ml/ |#> ? (1) 

This question concerns the nature and the validity of the approximations that 
deliver eq. ([[]) starting from the solutions of the Wheeler-DeWitt equation. There 
are two different (but related) approximations involved in this reduction process. 
First, the gravitational wave functions must be described by WKB solutions so as to 
have their phases governed by the gravitational Hamilton- Jacobi action. Secondly, 
the spread in energy- momentum of the state should be small enough so as to 
legitimate the description of its evolution by the Schrodinger equation 

idt\*) = Jd\N»H™{g l3 )\m) (2) 

in the gravitational background solution of eq. ([[]). The 4-D metric has been 
split in the 3+1 decomposition Q]: iV M is the lapse-shift 4- vector and g^ the metric 
of the spacelike hypersurfaces. T~C™{gij) = y/g~T® is the energy- momentum density 
operator acting on matter states. 

To address the validity of this reduction to a single geometry, one should in prin- 
ciple first construct matrix elements of in an enlarged framework in which each 
matter state is entangled with its own gravitational wave. Indeed, the Wheeler- 
DeWitt equation implies that each matter state determines its own gravitational 
wave (up to the specification of the quantum state of the linearized gravitons state 
that might be considered as part of the matter states in the present discussion). 
Only then one can ask in which circumstances it is legitimate to replace these 
entangled waves by a single one valid for all matter states. 

At this point it should be stressed that this procedure is not the one which 
has been generally adopted in the literature||-[T2|]. Indeed, in the "conventional" 
treatment, it is a priori assumed that the solutions of the WDW equation can be 
factorized into a single WKB wave function describing the gravitational background 
and something which is identified, to the lowest order in the inverse Planck mass, as 
the matter wave function evolving in this background according to the Schrodinger 
equation^. Higher order terms determine the gravitational corrections to the matter 
propagation ]TU|, |T2"fl . The main problem of this approach is that a part of these 
corrections is of purely classical character and directly follows from the choice of 

1 See in particular the recent article : "In this way the complete solution of the Wheeler- 
DeWitt equation is split into a product of a purely gravitational piece and a mixed piece repre- 
senting quantum field theory on that background." 
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the gravitational wave which has been factorized. In view of the present difficulties 
in the interpretation of the solutions of the WDW equation^, [13], [nf , one should 
carefully determine the origin of the corrections: Are they intrinsic to the WDW 
equation or are they induced by the approximations that have been applied to 
its solutions? This question is particularly relevant upon considering violations of 
unitarity@01|. 

To clarify these issues, another procedure has been adopted in |T5], |T6|, |T7| . The 
solutions of the WDW equation have been expressed as entangled superpositions 
of matter and gravitational waves and the evolution of the coefficients of this de- 
composition has been analyzed without factorizing a gravitational wave common 
to all matter states. After having obtained the dynamical equation governing this 
evolution, a gravitational background can be introduced by performing a first order 
expansion in the matter energy change. Then, the new expression of the dynamical 
equation reduces to the Schrodinger equation. Moreover, in order to minimize the 
corrections with respect to the exact background-free evolution, the background 
must be driven by the energy of the matter states under investigation. Finally, no 
violation of unitarity is produced by the passage to the new description based on 
a background. 

In the present paper, we pursue this analysis which was restricted to mini- 
superspace. Our aim is to explicitize the specific aspects introduced by the im- 
plementation of a Background Field Approximation (BFA) to gravity. In order 
to separate these aspects from the problems associated with the quantization of 
gravity, we shall work in a classical framework. However, we shall lead our inves- 
tigation so as to shed light on the quantum problem which consists in obtaining 
eq. ([TJ) starting from the solutions of the WDW equation. To this end, we shall 
analyze its classical counterpart: Starting from the Hamilton- Jacobi action which 
describes the entangled evolution of matter and gravity, what is the nature and the 
validity of the approximations that lead to the matter action satisfying the usual 
Hamilton- Jacobi equation in a given background whose evolution is fixed? Specifi- 
cally, we shall phrase this question in path integral terms by comparing neighboring 
Hamilton- Jacobi actions at fixed gravitational end points but differing in matter 
energy content. Having chosen this comparison, to apply a BFA to gravity is a well 
defined mathematical procedure. 

To first order in the change of the 3-geometries interpolating from the initial to 
the final condition, the linear change of the gravitational action equals the matter 
energy change in the background geometry. This is how one recuperates the usual 
HJ action governing matter propagation in a given background starting from the 
full HJ action, solution of the constraint equations Jip^avity _|_ _ g. rj^g ma ^ er 
energy term / d 3 xN tM H.™ (c.f. eq. (0)) is delivered by the linear "recoil" of the 
gravitational / pqdt term. We recall that the full HJ action contains only the sum 
of the / pqdt terms for gravity and matter. This recovery of the matter (light) HJ 
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action is not specific to gravity+matter systems. It occurs whenever one applies 
a BFA to the heavy sector of an enlarged entangled system composed of "heavy" 
and "light" degrees of freedom, see fl8|| . 

This classical result is just what we need in quantum cosmology. Indeed, 
the matter part of the full kernel, solution of the WDW equation, will obey a 
Schrodinger equation when the gravitational waves are WKB and when a first 
order expansion in the gravitational change makes sense. 

The validity of this first order expansion is controlled by higher order terms. 
These determine the gravitational corrections to the usual HJ matter action and 
their importance establishes a posteriori whether it was legitimate to treat gravita- 
tional degrees of freedom as "heavy" . They may be analyzed in classical settings. 
This crucial point deserves more explanations: The second order term has a struc- 
ture of the form Ae 2 x d^Sc where Sq is the gravitational HJ action evaluated 
on the background. However, in quantum cosmology, the origin of the spread in 
energy Ae may be purely quantum. Therefore, the structure of these corrections 
is fully determined by the classical analysis whereas its normalization may not be. 
This explains why we shall "recover" expressions that have been obtained in the 
quantum analysis^], [Tj|. Moreover we shall show that the corrections to the WKB 
gravitational wave functions are governed by expressions similar to those governing 
the BFA. This legitimizes the implementation of a BFA to WKB variables. 

Finally, the analysis of this classical problem has its own interest. It makes 
bridges between the techniques used when studying growth of local fluctuations 
(in particular as phrased by Steward and Salopek|B|) and those used in quan- 



tum cosmology. Moreover it also sheds light on the "issue of time" in (quantum) 
cosmology |T3|, pp| . Our analysis indeed does not coincides with the similar investi- 
gation performed by Barbour |21[ . 



In this article, we proceed in three steps. In Section 2, we analyze the action in 
the simple case of an homogeneous matter field in minisuperspace. Furthermore, 
we restrict ourselves to matter hamiltonians such that the amplitude of the field 
stays constant. Then the implementation of the BFA can explicitly and easily be 
performed in gauge independent terms. 

In Section 3, we work in a minisuperspace model in which the background is 
characterized by many degrees of freedom and we no longer put restrictions on the 
matter hamiltonian. Moreover we make use of techniques very similar to those 
required by the analysis of the solutions of the WDW equation. 

In Section 4, we generalize these results to arbitrary 3-geometries. We shall also 
comment on the treatment of second order corrections to the BFA which has been 



recently presented in (12 
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2 The background field approximation when mat- 
ter is characterized by a constant of motion 

As a warming up, we first show how a background field contribution is chosen and 
extracted from the HJ action describing the entangled evolution of matter and 
gravity. Then, we show how this choice determines the nature of the corrections 
to the description of matter evolution. To explicitize these corrections, we begin 
this Section by the analysis of the HJ action of matter and gravity without making 
any approximation, i.e. in a background free description. In order to simplify this 
analysis, we shall limit ourselves to cases in which the matter energy is determined 
by a constant of motion. 

We recall that in this Section and the next one, we work in minisuperspace. 
For these geometries, the metric element can be written as 

ds 2 = -N 2 (Od£ 2 + a 2 (Od 2 n 3 (3) 

where iV(£) is the lapse function, a(£) the scale factor and d 2 Q 3 the constant line 
element of the homogeneous three surfaces. 

The simplest matter system whose evolution is characterized by a constant of 
motion consists of an homogeneous distribution of comoving oscillators. Classically, 
these oscillators can be represented be the zero-component momentum of a massive 
scalar field [[L5[]. Their action is 

= J d^{p^-N{i)(pl/2 + M 2 4> 2 /2)} (4) 

where M is the frequence of the oscillations. On the equation of motion, their 
energy with respect to £ is 

N{i)e = N(£)M\A M \ 2 (5) 

where A M is the constant amplitude of the field. 

Another example consists of conformally coupled scalar photons of fixed (confor- 
mal) momentum k. On the equation of motion, their energy is Ne 1 = N^A^ 2 /a, 
see ||15|| . These examples can be generalized by considering non quadratic hamil- 
tonians in the field amplitude or by considering oscillations with an arbitrary a- 
dependent frequency u(a) treated in the adiabatic approximation so as to guarantee 
that their energy is iV(£)e m (a) = N(£)uj(a)\A m \ 2 . One can of course also consider 
a sum of these systems. 
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The background-free description 

The residual Bianchi indentity (see Chap. 27 in |3]]) requires that matter satisfies 

d a e m (a) = -47rP(a)a 2 (6) 

where P is the pressure. In all the above examples, this equation is satisfied because 
e m (a) characterizes on-shell matter propagation. Of crucial importance is the fact 
that the matter energy e m (a(£)) depends on £ through a only, thereby allowing eq. 
(@) to have this intrinsic lapse-free writing. 

In minisuperspace, the gravitational hamiltonian is given by 

N(OHg = {-G 2 *l + ™ 2 + Aa 4 ) (7) 

where G is Newton's constant, 7c a the momentum of a, k is equal to 0, ± for flat, 
open or closed three surfaces and A is the cosmo logical constant. Thus when the 
matter energy e m (a) is given, the HJ constraint equation, H G + H m = 0, reads 

- G 2 (d a S G (a, e m )f + kcl 2 + Aa 4 + 2Ga e m (a) = (8) 

The lapse N has also disappeared from this equation. In this Section, it will not 
be re-introduced thereby guaranteeing a gauge independent description. 
The HJ action S G , solution of eq. (§), is 

ra,2 

5 , G(Q2,e m )=/ dair(a,e m ) (9) 
where it (a, e m ) is the on-shell momentum driven by e m (a) 

7r(a,e m ) = -G^^na 2 + Aa 4 + 2Gae m (a) (10) 

Notice that the minus sign characterizes expanding universes. It arises from the 
unusual sign of the kinetic gravitational energy in eq. (§). 

The total HJ action (gravity + matter), solution of the constraint equation, is 
thus 

ra-2 r4>2 

S T = S G (a 2 , e m ) + S m ((f> 2 , e m ) = dair(a,e m )+ d<p p^((j),e m ) (11) 

where p^{4>, e m ) is the matter momentum, determined by the energy condition 
H m {P(j>,4 ) : a ) — e m( ( 3'). For definiteness and simplicity, from now on we consider 
only the case of harmonic oscillators. In this case, one has P(/,(<f>, e) = y/2e — M' 2 <f)' 2 , 
see eqs. (^|, 0) and compare it with eq. ( jID|) . Conformal photons or more general 
oscillators are treated along similar lines. 
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At this point we have imposed that the total energy vanishes and that both 
gravitational and matter propagations occur along classical orbits. However, when 
working with fixed the end point conditions (aj., 4>i, a 2 , 2 ), St is not fully extremal: 
It should still be extremized with respect to variations of ef20|, |23| . This stationary 
condition determines the "saddle point" value e = e(ai, 4>i] a 2 , 2 ) and reads 

a.2 r<t>2 



ra 2 r<p 

— d e / da n(a, e)\ €= i = d e 

J ai J (pi 



a 2 n r<f>2 1 

da 7r^ = / W—u^ (12) 

ai &Tc(a, e) Jfa p<t>{(p,€) 
t(a 2 ,e) = t ces i urn ((p2,£) 

On the left hand side, the gravitational integral has defined t(a 2 , e), the proper time 
lapse from a\ calculated from the propagation of a driven by the "saddle point" 
energy e. 

On the r.h.s. the matter integral has defined the "cesium" time t ces i um ((f)2,e). 
In the case of massive harmonic oscillators at rest, it equals to the period 2tt/M 
times the number^ of periods. (In the case of conformal photons, one would have 
obtained respectively the conformal time on the l.h.s. and the number of conformal 
periods on the r.h.s., c.f. P^fl .) 

Is it through the equality of these times that one obtains the orbit 2 = 0(a 2 , e) 
characterized by e and which passes by at g^. Once this correlation is obtained, 
it can be recast in the usual way parametrized by t: 0g(t) , ai(t) . In this writing, 
the dependence in e is "blamed" on the dynamical variables. 

The unusual aspect in this derivation of the trajectory is that t has arisen a 
posteriori through the dynamics of both matter and gravity. This is because the 
usual linear term in Et exte rnai was not present in the total action St, eq. ([TT|) . 
Rather both Sq and S m depend non-linearly on e. Up to now indeed Sq and S m 
have been treated on the same footings. By the implementation of a BFA, this 
symmetry will be broken since the BFA will concern only Sq- In other words, only 
a will be treated as an "heavy" degree of freedom. Then, the notion of an "external" 
time to matter (light) dynamics will be recovered from the heavy character of a. 
Notice that the division into light and heavy variables is not a question of taste: 
The physical circumstances decide whether it makes mathematical sense. Notice 
also that this procedure based on a division between light and heavy variables does 
not coincide with Barbour's approach ||21|| to time in cosmology. In his approach, 



time is a redundant variables which depends on all degrees of freedom. 



2 There might be a discrete number of solutions to eq. ( |l2|) since we are dealing with periodic 
matter systems. However this multiplicity does not affect what follows since we can always choose 
to work with the lowest value. 
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The background field approximation 

We now apply a BFA to the gravitational part of the full HJ action Sq + S m . 
To this end, one should first choose the reference background. To choose it, one 
needs to focus on a given "physical" problem. Indeed the choice of the background 
should be tided up to the problem under investigation and made so as to minimize 
the errors induced by the approximation process. This procedure has been applied 
in quantum cosmology |l(|: Upon computing a given transition amplitude from a 



initial matter state to a final one, it was shown that the "best" background is 
determined symmetrically by the energy contents of the initial and final states. 
Any other choice leads to greater systematic deviations with respect to the exact 
(background free) description. 

We use the same philosophy in the classical framework: We want to determine 
the simplified description of matter propagation in a single background geometry 
when the final values of cf) 2 are centered around a given (p 2 and when the three other 
conditions {a\, <f>i; a 2 ) are held fixed. This imposes to work with the background 
determined by the end point values (ai, <f>\] a 2 , <fi 2 ) since any other choice leads to 
systematic deviations. To obtain the simplified matter description, we expand only 
the gravitational part of the HJ action in power of Ae = e— e B = e— e(ai, <j)\] a 2 , <p 2 ). 
(From now on, we use the subscript B to identify quantities evaluated on their 
background value.) Notice that we have kept e free since the new description of 
matter evolution still involves the extremization of St with respect to it. 

By developing S G to first order in Ae = e — e B , one obtains 

£0(02, e) +S M {<t>2,e) = S G (a 2 , e B ) + / da (e - e B )d e 7r(a, e)\ e=E B + S M ((f> 2 , e) (13) 

J ai 

As in the absence of approximation, Sg + S m must be stationary with respect to 
changes in e. Then, as before, d e Sx = determines the saddle point energy e. 
However, the new expression of this condition now reads 

t(a 2 ,e B ) = l^dcj)— l — (14) 
Hi P4,{<p, e) 

in the place of eq. (|T2D. The only difference with that equation is that e on the 
l.h.s. has been replaced by e B for all values of 4> 2 . Thus t(a 2 ,e B ), the background 
time, is external to matter dynamics and the orbit directly expressed as (f> 2 = <pi{t). 

The errors induced by this BFA are of two kinds: First e, solution of eq. (|I~4"D, 
differs (slightly when a is heavy) from the solution of eq. (|T^) and secondly, t(a 2 , e B ) 
is not the time in the actual geometry. This is the price to pay when one has chosen 
to work with a given background common for different values of <p 2 . It should be 
noticed that neither error affects the orbit (j) 2 = 4>e{t) since both are concerned with 
the choice of the orbit through e and the parametrization of a 2 by t(a 2 ,e B ). This 
absence of deformations of the orbit follows from the existence of the constant e 
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which allows to separate the action St into two unconnected pieces. This will no 
longer be the case in the next Sections. 

Eq. (|T4| ) has been obtained from the HJ action, eq. (|TT|), in three steps. One 
has first chosen the reference background driven by e B . One has then expand 
the gravitational part of the action to first order in e — e B and finally requires 
extremization with respect to e. This is how the stationary condition^ determining 
matter propagation in the background a B (t) = %s(t) is obtained from the total 
action Sq + Sm governing the entangled dynamics of matter and gravity. 

This procedure generalizes the usual test-particle (or test-field) approximation 
for two reasons. First, in that approximation, there is no influence of the matter 
energy on the determination of the background evolution. Here instead, the "mean" 
matter energy e B does determine a B (t), see eq. ([10]). Secondly, because of this 



partial backreaction effect, the corrections to the BFA differ from those to the 
test-particle approximation, see below. 

The validity of the background field approximation 

The validity of the approximate description based on eq. (|T4l) is governed by the 
influence of the higher order terms in Ae. The quadratic term of the expansion of 
Sg{o-2,^) around the background contribution Sc(ci2,£ B ) is 

(Ae) 2 (Ae)2 (Ae) 2 p daa dM^e) 

— — d e S G {a 2 ,e)\ e= - € B = ——d e t{a 2 ,e)\ e= - e B = —— — -| e=e - s 

I 11 J a, Cr 7r z (a, e 



ft B (a 2 ) (I 



GA 



e 



Ka B {t)+A[a B {t)} 3 + 2Gi 



■B 



(15) 



The second expression makes clear that this correction is governed by the "sus- 
ceptibility" of the background geometry. It is indeed controlled by the amount 
of change in the background time lapse dt = ded t t induced by the change in the 
matter energy de. 

In the last expression, we have factorized the contribution of the linear term 
(dtAe) so as to identify the relative correction to it. In order to give an estimate 
to this correction factor, it is appropriate to consider matter dominated universes 
for which Ge B >>> na + Aa 3 . Then the correction factor is simply Ae/4e s : the 
change in the matter energy divided by the matter energy e B which drives the 
background. Notice that this ratio is independent of G. Thus the BFA is not an 
expansion in powers of G. Similarly, in quantum cosmology, we shall see that the 
BFA is not an expansion in the inverse Planck mass. On the contrary, the test- 
particle approximation can be viewed as the linearized theory in (Ge) since e does 
not enter into the determination of the background. Eq. ( [T5| ) also shows that the 

3 In the language of MTWj|, it is called the "condition of constructive interference" in ref- 
erence to the machinery of quantum mechanics. Upon working in quantum cosmology, we shall 
see indeed that both eqs. (|l^, appear with this status. 
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best background geometry to describe matter propagation when the values of </> 2 
are centered around </>f is determined by e B = e(a 2 , <fif). 

To determine the influence of the quadratic term, one should add it to the 
l.h.s. of eq. (|i4|) and compute the new value of the saddle energy e. One finds 
that its effect is to shift the saddle point energy e by (e — e B )d e t B / d e t cesium . Thus, 
when is a light and rapidly changing variable when compared to a, this shift 
is negligible since d t t B <K d e t cesium . This inequality shows that the value of the 
correction term does not control alone the validity of the BFA for describing matter 
evolution. Indeed this validity is governed by a ratio in which the lightness of (f) 
also intervenes. 

The fact that the change of e is controlled by a ratio of "specific heats" is 
reminiscent to statistical mechanics: Consider a microcanonical ensemble of energy 
E composed of two systems 1 and 2 with densities of states given by Vt{E\) = e s ^ El ^ 
and uj(E2) = e s ^ E2 \ When evaluating the density of states of the total system given 
by / defl(E — e)u(e) by a saddle point approximation, one obtains the equality of 
the temperatures: d e S(E — e) + d e s(e) = which is the equivalent of eq. (|T2"D, see 
App. B in [|TjJ for more details. Moreover, close to equilibrium and when E — e^> e, 
one can treat the big system 1 at the BFA, i.e. develop S(E) to first order in Ae. 
Then the errors induced by this canonical approximation are governed by the ratio 
of the specific heats d^S/d^s. For homogeneous systems this ratio scales like the 
ratio of the volumes. Therefore the BFA can be considered as a large reservoir 
limit. 

We conclude this subsection by mentioning an interesting question [2"1|: What 



happens close to a turning point? In this case the quadratic correction term, 
eq. (|T5D, becomes unbounded since the gravitational momentum n(a, e B ) vanishes. 
This signals the "instability" of the parametrization of matter evolution by t(a, e B ). 
Thus, close to a turning point, t(a, e B ) is not a good parameter. A very good 
parameter is t(ir a , e B ) obtained by performing a Legendre transformation to Sq so 
as to fix the end value of the momentum rather than a 2 . A similar situation occurs 
upon considering the validity of the WKB approximation: Close to a turning point, 
this approximation becomes also completely unappropriate. However the physics 
is not directly affected since, it is the state (the ket) and not its expression in the 
position representation which governs physical processes through matrix elements. 

Application to quantum cosmology 

The application of these results to quantum cosmology is straightforward if one 
postulates that gravitational waves are WKB. In this case indeed, their phases and 
norms are determined by the gravitational HJ action, eq. (H). At first sight it 
might seem that the validity of the WKB approximation introduce new restric- 



tions. However, as pointed out in M and further clarified in |15j, the conditions 



legitimizing the use of WKB gravitational waves are concomitant with those gov- 
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erning the validity of the BFA: Both approximations are valid when the relative 
change of the kinetical energy of the scale factor is negligible. 

Upon quantizing matter systems governed classically by constants of motion, 
one obtains stationary eigenstates labeled by a discrete parameter n {e.g. the oc- 
cupation number) and characterized by an eigenvalue: 



H m (a)\n) = e(n)\n) 



(16) 



In the time-free description, stationarity means that d a \n) = 0. This implies that 
the kernel, solution of the WDW equation 



(H G + H m )K,(a, <p\ ai 







(17) 



can be decomposed as 



K(a, 4>; a l7 0i) = lC n (a; ai) K % 



n{<P, 91 



wherein the matter kernel K n (<j); (f>i) is equal to (<f)\n) (n\<pi) as in Schrodinger set- 
tings. Moreover since the structure of K, is still a sum of products, the BFA will 
concern only the gravitational kernels K, n (a;ai). This is true only for the simple 
mini-superspace cases we are dealing with. 

The gravitational kernel K. n (a; ax) is a solution the quantized version of eq. (^j), 

i.e. 

(G 2 d 2 a + na 2 + Aa 4 + 2Ga e(n)) /C n (a; a x ) = (19) 
In the WKB approximation, for a > ai, it is equal to 



K, n (a; ai) 



cxp 



ij^ da' 7r(a',e(n)) 



^n(a, e(n))7i(ai, e(n)) 



(20) 



We now determine the behavior of kernel K.(a, <p; a 1? <pi) when a is equal to 02 
and close to 0^. To make contact with classical mechanics, we use a saddle 
point approximation (i.e. a constructive interference condition) to estimate the 
summation over n in eq. (Jl8|). When using the WKB expressions eq. (0), the 
equation which fixes the saddle point value n B is equal to that one would have 
obtained by treating gravity at BFA if the background was the solution driven 



by the saddle point energy e 



B 



e(n 



This is clearly exhibited by factorizing 



fcn B {Q'2] O'l) and developing the gravitational kernels in An 



n 



n B 



K,(a 2 ,4> B ]a 1 ,4> 1 ) 



exp 



7r(q 2 ,e(n i? ))7r(ai,e(n i3 )) 
^ 7r(a 2 ,e(n))7r(a 1 ,e(ri)) 

da 1 7r(a, e(n)) — 7r(a, e(n B )) \ 
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IC nB (a 2 ;a 1 ) x ^ K n((f>2 '; <h) \ 1 + 



exp 



-2 



7r 2 (a, e(n s )) 

/;- M * W «)-^)}{i+o(^)} 



(21) 



where A n e = e(n) — e(n ). The stationary phase conditions are respectively 



"2 



<9 n ^/ davr(a,e(n))+^ In AT n (0f ;<f) 1 )\ = (22) 

ai 



de{n) 
dn 



t nB (a 2 )+td n \nK n (^; ( j )l ) = (23) 



Eq. (p2|) is obtained from the first expression of eq. ( pi]) and does no rely on 
any BFA. It can be viewed as the quantum version of eq. ([12]). Instead eq. (p3|) 
follows from the second expression in which a BFA has been already performed, 
c.f. eq. ([14]). This later condition is the usual "constructive interference condition" 
obtained from QFT in the background characterized by t n B(a). By construction, 
these conditions coincide since the expansion in An was performed precisely around 
the saddle point value n B . 



From this agreement and the former classical analysis leading to eq. (|14D, one 
immediately deduces that the saddle point phase of the kernel considered as a 
function of 02 in a neighborhood of <p B is equal to the BFA phase when non linear 
terms in n — n are discarded. (Notice that these correction terms are governed by 
A n e/7r 2 as in the classical case.) Moreover, when they are discarded, the spread in 
n — n agrees with the BFA spread since in both cases (with or without BFA) it is 
determined by the (second derivative of the log of the) matter kernel K n ((j); 0i) only. 
Notice that this spread is intrinsic, i.e. determined by the physical circumstances 
and h. This is to be opposed to the classical spread Ae in the former subsection 
which was fixed by the "external" choice of 02- However, the mathematics involved 
are the same. 

The above three agreements imply that /C(a 2 , 4>2] ai, 4>i) satisfies the Schrodinger 
equation in the background geometry a B (t) for 2 sufficiently close to the back- 
ground value 4> B so as to legitimize the neglection of higher order terms in An. 
Notice however the unconventional property of this Schrodinger equation. Since 
the background geometry is dynamically determined by n, c.f. eq. (|2"2"|), it changes 
when a,2 varies except if one tunes <j) B so as to stay on the orbit <p2 = <P{(i2,n B ). 
Because of this, it is appropriate to work with fixed n rather than with fixed end 
point values on the matter field. This is what has been adopted in [[1^, [T7|] to study 



transition amplitudes in quantum cosmology. This is also what we shall adopt in 
the next Sections in the definition of the background. 

As in classical mechanics, higher order terms in An (or Ae) induces corrections 
to the determination of the saddle point value n. To determine their physical 
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consequences it is mandatory to introduce interactions leading to transitions among 
states, see [16 , [17]] • The main lesson of these works is that it is not necessary to 



first factorize the mean gravitational kernel and then search for the description of 
matter evolution|. Indeed, the WDW equation determines their evolution directly 
in terms of a, i.e. in a background free description, see eq. fllBp for a simple 
example. However, as in eq. (pip, one can always a posteriori perform a BFA (an 
expansion in n — n B ) in order to obtain a more conventional description formulated 
in a single background |[20||. Doing so one induces non linear terms in n — n B whose 



role among other things is to introduce phase shifts governed by d^Sc- 

In resume, eq. fl2T|) shows that the corrections to the kernel with respect to its 
Schrodinger BFA expression are of two kinds: those due to quadratic and higher 
order terms in An and those due to the WKB approximation, eq. (|20|). In a next 
paper, we shall analyze the consequences of both types of corrections on possible 
violations of unitarity. 



3 The background field approximation for gen- 
eral mini-superspace models 

In this Section we no longer impose that matter dynamics are governed by a con- 
stant of motion. Therefore we are obliged to apply the BFA to extremized actions 
since we have lost the possibility of applying a BFA before having extremized the 
total action with respect to e, see eq. (filf ). There are two ways to implement this 
approximation. 

The first one consists in keeping the lapse function N(£) and working with 
trajectories parametrized by £. The main advantage of this approach is that the 
whole space-time structure of the background emerges upon applying the BFA to 
gravity. Indeed since both a B (£) and N B (£) are determined, one obtains matter 
propagation in a curved background which possesses a given time parametrization. 
A second advantage is that it offers the possibility of considering "off-shell" matter 
configurations. At first sight it might seem completely absurd to consider the 
gravitational reaction to off-shell matter configurations since r.h.s of Einstein's 
equation will be incompatible with the Bianchi identity, i.e. eq. (0) will not 
be satisfied. However, in the path integral, most of the matter configurations 
from 0i to 02 do violate this equation. Therefore it is a legitimate question to 

4 It is even a nuisance! The reason is that transition amplitudes are independent of the mean 
energy: they only depend of the energy of the two states characterizing the transition. However, 
when searching for the description of these transitions in the background geometry driven by the 
mean energy, one obtains series in which all terms depend parametrically on this energy. This 
is nothing but the consequence of having first factorized the mean gravitational kernel. By no 
means it is an intrinsic properties of the solutions of the WDW equation. 
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ask whether the action governing these configurations can be conceived as arising 
from an enlarged action governing the evolution of both gravity and matter. This 
analysis is performed in Appendix B. 

The second approach is more intrinsic since it is based on the HJ action rather 
than on trajectories. It does not make use of the lapse function and furthermore 
proceeds along lines very similar to those used in quantum cosmology. Therefore, 
we shall present here this second approach. The first one is developed in Appendix 
A. The reader unfamiliar with HJ techniques is invited to first read this Appendix 
which proceeds along more conventional mechanical lines. 

In order to prepare the analysis of the next Section which deals with gen- 
eral 3-surfaces, we work with backgrounds described by many degrees of freedom, 
X 1 . These consist on the scale factor a, on the homogeneous part of some mas- 
sive inflaton field and possibly also on gravitational Fourrier modes highly ex- 
cited. The multidimensional character of the background trajectories will introduce 
new difficulties since there are now "transversal" directions to the background or- 

bitsQ[0,0,[g. 

The total hamiltonian we are considering is of the form 

N(OH T = iV(£) [HxfaX*) +H m (4>,p^,X i )] (24) 

We only impose that H x and H m be quadratic in the momenta 7Tj and p^. Notice 
also that we have not included in H m a "non- minimal" coupling to 7T; . The inclusion 
of this additional coupling presents no difficulty however. 
The total action satisfies the HJ constraint equation 

^{X k )d x ,S T d xi S T + V(X k ) + tf m (0,<V? T ,X l ) = (25) 

As in eq. (|TTD , St can be written as a sum of terms 

■ ( IX'n,iX'.{Xi.n 2 \) + I 



S T = S X + S m = J* 2 rfXV^r, {X 3 2 , 2 }) + f 2 dfotih {X J 2 , 2 }) (26) 



where each on-shell momentum is parametrized by its conjugate position. In brack- 
ets we have indicated that they depend parametrically of the end point values of 
all variables. It is on this dependence that the BFA will act. 



5 It might seem a little bit antinomic to consider a plurality of backgrounds. Moreover, there is 
a simple way to reject this plurality: It consists in treating all but one coordinates X % (say e.g. a) 
as "matter" degrees of freedom. One would then recover a situation similar to that of the former 
Section in which the "matter" energy univocally specifies the orbits of the remaining background 
coordinate. This is effectively what has been adopted by Halliwel and Hawking in ||. In what 
follows however, we shall proceed with background orbits characterized by many coordinates X 1 
in order to confront the new difficulties. 
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Before proceeding to the implementation of the BFA, we should carefully de- 
termine the background(s). In the former Section, this was easily achieved by 
the specification of the constant of motion e B . As we discussed in the "quantum 
cosmology" subsection, this amounted to work with the end value of evaluated 
along the orbit characterized by e B . In the present situation without constant of 
motion, we define the background by fixing the initial value of the matter energy 
H m (ai,</>i,p^), i.e. we fix the initial momentum p^i = p B . Then the end value of 
the 4> field is evaluated along the orbit specified by p B as well as the end-point val- 
ues X\. Indeed the specification of the matter energy no longer suffices to specify 
univocally the orbit. Therefore the set of X\ will both determine the time lapse 
along the trajectory (as a did in the former Section) but also determine the orbit, 
i.e. the location in the "transversal" directions. 

With these aspects in mind, we now search for the approximate description of 
matter evolution when the change in 5^ = iii(X'\ {X 2 , 2 j) — ir^X 1 , {X 2 , p B j) 
induced by the change in the matter configurations is treated to first order. To 
this end, we write 

S T (Xi 2 ) = S B (Xi) + AS(Xl, 2 ) (27) 

where S^X^) is the background action given by f x i dX l ix i (X' 1 , jx 2 , p^j). By 
inserting eq. (|27j ) into the constraint eq. (^5|) one obtains 

C^(X k )d x ,S B d xi AS + H m (<j>, d^AS, X*) - 

+ ho ij {X k )d xi AS d X3 AS = (28) 

where we have used the fact that S B is the solution of eq. (^5]) driven by the 
background matter energy H^X*). Notice that H^X*) is orbit dependent but 
nevertheless acts as a potential energy along each orbit. Moreover, along each 
orbit, the proper time lapse to go from X{, <pi,p B to X\ satisfies (by definition) 

C ij (X k )d x *S x ' d x d B {X l ) = 1 (29) 

This is just another way to express the relationship between velocities and momenta 
which is usually written as dX lB /dt = C^irf when one has t at our disposal from 
the outset. Here instead, t B (X l ) is re-introduced and defined by —OeS x where dE 
is an infinitesimal constant change of the matter energy H^. Then by varying eq. 
(|25| ) applied to S x with respect to dH B = dE, one obtains eq. (p9|). 

In the present formalism, to apply a BFA consists in neglecting the last term 
of eq. (PS]). This amounts to assume that d X iAS -C d X iS B , i.e. 57r» irf. Then, 
using eq. (p9|), the first line of eq. (p8|) reads 

d t AS^,X lB (t)) + H m ((p, d^AS, X iB {t)) - H B (t) = (30) 
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Thus, to first order in 57Tj, AS satisfies the HJ equation for the matter field in the 
background determined by X\,p B . Then, as usual, it can be decomposed as 

AS(M = J rf# (0,{0 2 }) - j dt'{H m (<P,P^X iB (t'))-H*(t')} (31) 



Together with eqs. ([26"1 , |2"7|) , this expression shows that the energy term —AH m dt 
originates from the linear recoil of gravity, 5iTidX l , induced by the change in the 
matter trajectory. 

Notice that all heavy degrees of freedom X % contribute to the determination of 



the background time. In this we agree with what has been advocated in [fZI] safe 
for the fact that the light variables only intervene in the definition of time through 
their "mean" energy H B . To our opinion however, the important lesson of eq. (p0|) 
is not that t B (X l ) is a "redundant" parameter but that it is precisely the one which 
does the job: to parametrize the dynamical evolution of the light variables that 
did not directly participate in its determination. The reasons for which t B (X l ) 
is the right parameter are the following. First, the change of the light trajectory 
that we want to parametrize induces through the constraint equation a change in 
the energy available to the heavy variables. Secondly t B (X l ) is defined by —DeS 3 -, 
i.e. the conjugate to an infinitesimal change of that available energy. Then to first 
order in the finite matter energy change, these two reasons imply eq. (|30D . 

Notice also that in spite of the fact that the new matter configuration defines 
a new background orbit, the "transversal" part of the change does not show up 
to first order in since it is projected out by the contraction with the gradient 
C l WiS B dj which is perpendicular to the surfaces S B = Const. This is a well known 
feature of first order perturbation theory in classical mechanics, see e.g. | f2~5f . 

We now verify that the neglected term in eq. (|30D governs back-reaction effects, 
i.e. the modifications of the matter dynamics induced by its own energy through 
the "susceptibility" of the background. To estimate this term perturbatively we 
use eq. ( |3"T| ) and we obtainP] 



^(X k )d x ,AS d X jAS = ^(X k )d x d B d X3 t B (AH m ) 2 (32) 



6 As pointed out to me by Serge Massar, eq. (|32| ) is exact only in the limit of slowly varying 
AH m (X l ), i.e. when AH m ^> (X^ — X\)dx^AH m . This is the usual condition appearing in 
adiabatic treatments [p2[ . Thus, for rapidly varying AH m or extremely long trajectories, one 
should exactly take into account the transversal dependence of AH m (X l ) associated with the 
change in the end point value X\. This amounts to compute the Jacobi fields of S^- The reader 
interested by the exact calculation will consult Section 5 of fjl| . Let us mention here that it is 
necessary to perform the exact calculation only when the variables X % are not sufficiently "heavy" . 
Notice also that our estimate, eq. (^2|), does not coincide with that of fU^ . In that article, the 
longitudinal is defined by the surfaces S% = Const, whereas here it is defined by t B = Const. A 
detailed comparison of these estimates goes beyond the scope of the present paper. For further 
discussion see next Section. 
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This term generalizes what we obtained in eq. (|i~5|). Indeed for backgrounds 
characterized by a single variable, one has C(X)dxt B dxt B = l/C(X)n x which is 
the integrand of the fourth expression in eq. (jT5|). Moreover, in multi dimensional 
cases, when the change in energy AH m is constant, eq. ( |3~2D still contributes to AS" 
given in eq. (^T|) by the addition of (AH m ) 2 d E t B /2. By d E t B we designate d E S B . 
To evaluate it, one should vary eq. (f29|) with respect to dH B = dE. One finds, 



(33) 



d E {c^{x k )d xl s B dxifix*) = 1} 
c ij (x k )d x d B d x d B - c ij (x k )d x ,s B d XJ d E t B = 

Using once more eq. fl29|) to rewrite the second term of eq. (0), one gets 



d E t B = f dt'C ij (X k )d x d B d x d B 
Jo 



(34) 



From this equation one verifies that eq. (^) introduces a source term to eq. (p0|) 
which integrated over t gives the announced result in the cases of constant AH m . 

When the last term of eq. (p8|) is replaced by the r.h.s. of eq. (|32|), one 
obtains an equation for the matter propagation since all X 1 are evaluated on their 
background values and parametrized by t. Thus one can view the last three terms 
of eq. ( ^81) as defining the "effective" matter hamiltonian. This hamiltonian gives 
rise to unusual Euler-Lagrange equation since it is quartic in = d^AS. This 
results from the "elimination" of the degrees of freedom X 1 . A simple procedure to 
handle this hamiltonian consists in treating the correction term perturbatively, in 
a manner similar to what is done in a Hartree approximation and also in agreement 
to what we have learned in the former Section. Namely, the quadratic correction 
term only shifts the saddle point energy e by a term proportional to d 2 SG/d 2 S m . 

In this procedure, one first uses the unperturbed action eq. ([H]) to obtain the 
energy difference AH m (t) evaluated along the unperturbed equations of motion. 
Having obtained this function of t (a c-number in the quantum mechanical lan- 
guage) one linearizes eq. ( p8|) with respect to H m ((j), d^AS, X l (t)) — H B (t) so as to 
keep the quadraticity of the HJ equation in d^AS. One gets 



d t AS 



H m (<f>,dtAS,X\t))-HZ(t) 



1 + -AH m {t)C l 'dd B d 3 t B 



(35) 



By construction, the correction term is now expressed as a factor. Moreover, since 
it depends only on t, it can be absorbed in the l.h.s. by a redefinition of t. Thus 
its sole effect is to modify the time-parametrization of the background orbit X lB {t) 
and therefore the lapse of time between the end points X\ and X l 2 . Notice that 
this procedure can be refined by determining self-consistent ly AH m (t). 



Application to quantum cosmology 
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The aim of this brief subsection is to show how similar are the quantized versions 
of the equations we just obtained. 

As in the former Section, in order to extract the background contribution, we 
write the total kernel, solution of the WDW equation (the quantized version of eq. 
Q55D), as a product of the background gravitational propagator times the rest: 



K(Xl 2 ; XI 00 = K B {X\- X\) x AK(X l 2 , 2 ; X{, X ) (36) 
The background gravitational kernel is the solution of 



-C l \X k )d x ,d X i + V(X k ) + H^(X k )\ TC B (X l 2 - X[) = (37) 
Since the variables X % are heavy, it is legitimate to use the WKB approximation: 
K^ KB (Xi- X\) = ^Det.d^d^S* exp {iS%(X* 2 )} (38) 

where the normalization factor is given by the enlarged Van Vleck determinant 
which occurs when one works at fixed energy, see f26|| . By X 2 and X", we designate 
the sets of variables (X 2 ,dE) and (X(,dE). (Notice that we should not worry 
about the normal ordering ambiguities of the operator C % \X k )d Xl d X j since we are 
working in the WKB approximation.) 

Using the fact that the background kernel is the WKB solution of eq. (0), one 
immediately obtains the equation for AK 

C ij (X k )d X iS^ id xi + H m {<P, zfy, X 1 ) - H B (X' 1 )] AK(X\ 0; X{, 00 

' -C i3 (X k )d xl d xi ] JDet.d x ,d x ^S B AK(X\ 0; X{, 00 = (39) 



The parallelism between the classical version eq. ( 28) and this equation is 



manifest. In the quantum version, upon neglecting the second line of eq. (p9|), 
i.e to first order in small gradients d X iAK and d X iDet., AK is the usual matter 
kernel and it satisfies the Schrodinger equation in the background geometry driven 

by H'JXX 1 ). 

Higher order terms govern the corrections to the WKB and the BFA approx- 
imations. Their structure results from the a priori factorization of JC B 7KB . Thus 
one should be careful in interpreting these terms as providing intrinsic corrections 
to the WDW equation. It should be indeed recalled that for unidimensional mini- 
superspace models, the intrinsic corrections have been determined from the exact 
background-free equation governing matter evolution [|17|| which was obtained with- 
out factorizing a given gravitational kernel nor postulating that the WKB expres- 
sions constitute a good approximation. It was show that the intrinsic corrections 
do not reproduce those resulting from a a priori factorization. The challenging 
question is thus to extend this analysis to the present multidimensional case. We 
hope to be able to report on this question. 
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4 The background field approximation in cos- 
mology 

In this Section we indicate how the results of the former one can be generalize to 
arbitrary three-geometries whose metric is given by g^{x). The total hamiltonian 
which governs the entangled dynamics of g'^{x) and some local matter field <f)(x) 
has the following structure 

H T = Jd i xN^^x)Hl(^x) = Jd 3 xN»(^x) {n^g^+H^A,^)} 

(40) 

The four-vector field iV M (£,x) describes how successive 3-surfaces are pasted to 
form a 4-dimensional space-time, see ||]. 7Y^(£,x) is the gravitational energy- 
momentum density. It generalizes eq. (0), it is quadratic in the gradients d x kg^(x) 
and in the conjugate momenta ^ (x) but is a non- linear function of g^(x). The 
matter energy-momentum vector density H™(£, x) depends on 4>(x), d x k<p(x), p<f>(x) 
and g^(x). 

The main difference with Section 3 is that there are now four local equations of 
constraint in the place of eq. (^) which result from the four Einstein's equations 
On^St = "H^iC, x) = 0. In the HJ guise, they read 

l -C i ^\g k \x))d g , ] S T d gifj ,S T + V{g k \x)) + H^{x), d^S T , g^{x)) = (41) 

-2{d gi iS T ) lj +HT{(f>{x),d ((> ST,g ii (x)) = (42) 

For further details concerning the kinetical matrix C or the potential V see M. 
The important point is that these constraints do not decouple since the potential 
term V(g kl (x)) contains gradients connecting the fields at different points. Upon 
neglecting this term one obtains an action which is just a sum of decoupled actions 



evaluated at each point. As shown in [19|, one can then perform a perturbative 



treatment in those gradients. What follows can be conceived as a generalization 
of these works in that the small parameter which governs the expansion is the 
change of the matter hamiltonian itself. Therefore, our (less explicit) treatment 
also applies to arbitrary large gradients. 

In spite of these gradients, the total action from the initial configuration g\ J (x), (f>i(x) 
to the final one g 2 (x) , 4>2{x) which is a solution of the constraints can be decom- 
posed, as in eq. (|26|) , as 

S T = [f X) dg^(x) -„(//": {£, <M) + d<j>(x) P4 ,(<j>; {g k \ M) (43) 

where the local on-shell momenta are parametrized by their conjugate field. In this 
paper, we shall not discuss the difficulties that might prevent to actually compute 
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St when the initial and final configurations are given, i.e. to solve the "thick" 
sandwich problem. Rather we shall only use the structure of St and the fact that 



it is a solution of eqs. ( |41| , [42| ) since it is all what we need in our implementation 
of a BFA to gravity. 

To conform ourselves to the "standard procedure" JTU], [T2|, [IB] (later on we shall 
comment on its appropriate character), we treat all g^(x) in the same footing as 
"heavy" degrees of freedom. Thus, as in the former Section, we choose an initial 
matter configuration (p<pi(x), 4>i(x)) which determines both the background matter 



density H B and the background trajectories. Then, as in eq. (|27| ), we decompose 



St as S^gY) + AS , (0 2 , g^) and we perform an expansion in the changes 57iij(x) 
which result from the change in the matter field configuration induced by the 
replacement of the background value (p B (x) by the new field configuration specified 
by the final value 02 (x). By inserting this decomposition in eqs. f42j), we obtain 
respectively 

C^\g k \x))d g ,S B d giV AS + H (<P(x), d^AS, g V(x)) - H B (x, g*(x)) 

+\c i ^Xg k \x))d^AS d giV AS = (44) 

-2(d gij ASy j + (Hr(<P(x),d<f,S T ,g t3 ) -HT B (x,g l ^x))) = (45) 

In order to relate these local equations to a single HJ equation governing matter 
propagation in a given background metric, we must introduce by hand the lapse- 
shift vector field N^(x, t). These four functions then define univocally a single time 
parameter. It is given by t = —8eSx where dE is an infinitesimal constant change 
in the total matter energy / d 3 xN' 1 (x)T-l 7 ^ B (x). More precisely, this "redundant" 
parameter satisfies 

d 3 x{N°(x)C i,j ' ij (g kl (x))d giV Si d g ,+2N\x)^d^}t B (g l \x),N^x)) = 1 (46) 

This equation is obtained by varying eq. (|4T)| ) with respect to dE and by imposing 
Ht = 0. It allows to recover the usual relation between the momentum and the 
velocity: dg»(x)/dt = iV°C^V iy + 2N^ . 

Using this time parameter and eq. (flop , and upon neglecting the second line 
of eq. ([P]) which is quadratic in the small gradients d g ijAS, the first line of that 
equation becomes the HJ equation 

d t AS = J x) {^(^(f, i), %,,)AS, x)) - H; B (t, i)} (47) 

evaluated in the curved geometry solution of the Einstein's equation driven by the 
"mean" matter energy density l~L™ B {x). 

In brief, the novelty introduced by considering arbitrary 3-geometries relies on 
the necessity of introducing from the outset the four functions N^(x) in order to 
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rigidify the propagation in super-space from the initial 3-geometry to the next ones. 
We emphasize that the dynamical justification arises from the presence of gradients 
d x g % \ d x ((> in the potential terms of Ht- 

In what follows, we shall not present the quantum version of this derivation 
which would lead to the Schrodinger equation in a curved geometry. It is straight- 
forwardly obtained by generalizing eqs. (|36| - |39D . (The interested reader will find 
this quantum treatment in [12] which is however restricted to empty solutions of 
Einstein equations, i.e. to 7Y^L = 0.) Once more, the reader will verify that the 
structure of the quantum version is fully determined by the implementation of a 
BFA in classical settings. 

We shall neither present the computation of the quadratic corrections which 



arise through the second line of eq. (23). They can also be found in [12]. What 



we shall do instead is to discuss on physical grounds the appropriate character of 
applying, as we just did, a BFA to all gravitational degrees of freedom. 

We shall base our discussion on the fact that the quadratic corrections to AS 
are given by a term of the form J J T^ lu D flUfJ-u T^i where D is a gravitational Green 
function. These interactions are in strict analogy with the current-current interac- 
tions in electro-magnetism obtained by integrating over the photon field. However, 



as discussed in ||12|| , the Green function which is obtained by treating all gravita- 
tional degrees of freedom on the same footing is rather unconventional in that it 
vanishes both on the initial and final 3-surfaces, in contradistinction with the usual 
retarded and Feynman functions in classical and quantum settings. This uncon- 
ventional feature results from the fact that the g l i have been completely specified 
on both the initial and final 3-surfaces. 

This unpleasant feature reveals that the "heavy" character that has been a 
priori attributed to must be called into question. Indeed, so far, g^ has been 
treated like the W gauge boson in the derivation of the four fermions weak in- 
teractions model. In that case, it is the rest mass of the W which legitimizes its 
elimination for low energy processes since no on-shell W is ever produced. For 
gravity, the Planck mass does not act like the W mass and no mass threshold pre- 
vents the production of on-shell gravitons. For this reason, it is inappropriate to 
functionally integrate out the gravitons. Indeed, in quantum gravity, these gravi- 
tons are correlated in phase and amplitude to the quantized matter sources. This 
quantum entanglement prevents the possibility of describing the final gravitational 
configurations by c-number functions giving rise to a classical background action. 
More intuitively, on which physical basis could one justify an unsymmetrical treat- 
ment in which rare (i.e. not coherent) on-mass gravitons are considered "heavy" 
and thus described by WKB waves on which one applies a BFA and photons with 
the same frequency are "light" and thus fully quantized? 

In order to treat gravitons and photons on an equal footing, one must therefore 
abandon the idealized procedure in which all gravitational degrees of freedom are 
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treated on the same footing. By treating the light part of the gravitational degrees 
of freedom like matter degrees of freedom, one "recovers" what has been adopted 
in cosmology to study the growth of perturbations 0, see also for a discussion 
on the dynamical role of light and fluctuating gravitational degrees of freedom in 
explaining why the heavy coordinates are driven by mean values. 

In conclusion, it is physically unjustified to search for a treatment of the solu- 
tions of the WDW in which all gravitational degrees of freedom are treated on the 
same footing. However what is lost from a geometrical point of view is compensate 
by a dynamical justification: in each situation, it will be convenient and appro- 
priate to treat some gravitational degrees of freedom as heavy but it will also be 
mandatory to keep the other ones fully quantized. The same division also applies 
to matter degrees of freedom. The implementation of this program requires further 
work which will hopefully be presented somewhere else. 
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5 Appendix A. The background field approxi- 
mation based on trajectories. 

Given the hamiltonian eq. (|24|), the total "off-shell" action we are considering is 

S T = 1^ ^{nX'+p^-N^) [H x {^X l ) + Hmifap^a)]} (48) 

In order to have univocally specified expressions, we must choose a "representa- 
tive" for iV(£). Indeed, at fixed end point values X\,<f>i to X|,0 2 , the action 
St is invariant under reparametrizations of £ which leave the end point values 
£ = 0, 1 fixed. The reparametrizations which satisfy this restriction leave the in- 
tegral Jo d^N(^) = N invariant. Thus one can only (but always) choose as a 
representative the constant lapse iV(£) = N. This choice is the one usually made 
when describing the path integral of a relativistic particle. In that case, ds = Nd£ 
defines the Schwinger fifth-time. Adopting the same choice, our "representative" 
off-shell action is 

rN 

S T = dt^X'+p^-lHxi^X^+H^p^X*)}} (49) 

J 

We have used the symbol t defined by dt = Nd£, since in this "gauge" of constant 
lapse, it indeed corresponds to the proper time. In this gauge, the total hamiltonian 
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is a constant along the classical trajectory which relates X\, (pi to X 2 , 02 m a proper 
time N since it is not explicitly time dependent. 

Since physical propagation must occur with vanishing total energy, one must 
impose 8nSt = Hx + H m = 0. This extremization condition determines the 
"saddle value" N = N(X 2 , 2 ; X{, <pi) which, in our gauge, equals t(X 2 ,(f) 2 ), the 
propre time lapse from X\ to X\ given (pi and (f> 2 . d^Sx = plays thus the role of 
d e S T = in Section 2: Both fix the proper time lapse on the classical trajectory. 

Thus the fully extremized action is 

S T = S X + S m = f N dt {mX 1 + p^} (50) 

J 

In this action, X l {t) and 0(t) are evaluated along the equation on motion, they 
thus depend parametrically on N. 

We are now in position to implement a BFA. As in Section 2, we shall linearize 
only the gravitational part of the action in order to obtain the simplified description 
of the matter evolution for values of 02 centered around the background value 
$2 ■ The difference with that Section is that we shall consider the difference of 
extremized actions since no longer can work with e off shell. 

Thus, we want to expand S^A^, 02, ^(X 2 , 02)) around the background action 
St(X2, 0f , N B = N(X 2 , 0f)) to a first order change of the background coordinates 
SX z (t) and the lapse SN. We proceed in three steps. 

In the first step, we replace N(X\, 2 ) by its background value N B = N(X\, 0f ) 
in the first action. To first order in SN, the action is left unchanged since we are 
working with vanishing total energy and we are not modifying the end point values 
of and a. To consider only the first order change in 5N is a valid restriction 
because, when X l {t) are heavy, i.e. when second order change in 5X l {t) are neg- 
ligible with respect to first order change, the same in true for the changes due to 
5N. In other words when X l (i) are heavy, iV is also heavy. (This is exactly like in 
statistical mechanics. The change of the temperature induced by a fluctuation of 
the energy repartition scale like the inverse volume of the heavy system.) 

In a second step, we consider the first order change in 5X l {t) and SiTi(t). We 
have 

S x (Xl 02, N B ) = Sx{Xi, <P B , N B ) + J o dt fax* - 8X%) + 5XV^ (51) 

Since we are working with fixed end points X[,X 2 , the boundary term vanishes. 
Moreover since we are on-shell, 

57i t X* - 8X% = 6 xi>ni (H x + H m ) (52) 

where ^x 4 ,^ means the variation of the quantity induced by the sole variations 
of X % and 7Tj, i.e. with 0(t) and p^t) evaluated along their background classical 
values. 
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In the third step, we re-use the vanishing of the total energy for both classical 
trajectories. Thus 



^x i ,^ i (Hx + H m ) 




The second equality follows from the fact that H x is independent of and p^. The 
last one means that the variations with respect to and p^ should be taken to 
all orders and not to first order only in the case of 5X l and 5iii. We emphasize 
that both terms in this last expression are evaluated in the background geometry 



The 02 dependent part of this action is the usual HJ matter action evaluated in the 
time dependent background X tB (t). The simplest way to verify it is to consider 
first order variation of the end point values X 2 . At this point indeed, we cannot 
vary the time lapse N B since it is fixed by a 2 . The first order variation of St with 
respect to an arbitrary variation of the end points X 2 gives 



5 S T = 5X1 d n S T = 5X1 { 9i {Xi) + d A N B (JC 2 , 2 ) H m (d*S T , 0, X*)} (55) 



where the dependence in X 2 which is irrelevant for the matter dynamics has been 
put in the functions Using N B = t(Xl), we can now determine the change of 
St when the 5X 2 are evaluated along the classical trajectory, i.e. when they are of 
the form 5X 2 = C l ^djS x 5t. By definition of t one obtains 



It is thus through changes of X % evaluated along the background orbit that one 
recuperates the time dependent HJ equation for the matter degrees of freedom. 

6 Appendix B. The gravitational response to off- 
shell matter configurations 

In the former Appendix we saw that the linear change of the background action 
Sx induced by the change of the matter configurations is equal to / dtAH m eval- 
uated along the background trajectory X lB (t), see eq. ([54]). However it should 





+ j o dt{ P4> (t)j)(t) -H m (p4t)A(t),X w (t))}(54) 



5S T = 6t d t S T = 5t {H m (d^S T , 0, a B (t)) + 9l (t)d t X w (t)} 



(56) 
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be noticed that this result has been obtained by comparing neighboring on-shell 
matter configurations. It would therefore be very useful to extend this result to 
off-shell matter configurations in order to be able to use this relationship in a path 
integral formulation of quantum mechanical matter propagation which does involve 
off-shell configurations. 

Thus the problem one should analyze in classical settings is to determine the 
gravitational response to a given arbitrary matter configuration </>(£) , Pj,(£) ■ Strictly 
speaking, this cannot be achieved since one cannot find a stationary gravitational 
action whose source would be this particular matter configuration. Indeed, there 
is no reason that this configuration satisfies the conservation equation eq. (||). 
Therefore, the Bianchi identity tells us that no (on-shell) background could be 
driven by this source. 

However this does not prevent us to characterize the departure from stationarity 
when the configuration 0(£)>2ty(£) is close to a classical orbit 4> B {£,),p B (£)■ As we 
shall see, this departure is just what we need to recover the usual off-shell matter 
action in a given background that one uses in a path integral formulation. To prove 
this result, we shall simply adapt the procedure we used in the former Appendix 
so as to take into account the off-shell character of the configuration </>(£), p^O- 

The total unextremized action is given by eq. (|49|) since one can still work 
with constant lapses N as representatives. (The off-shellness of </>(£) does not 
spoil reparametrization invariance). Therefore, at fixed N, one can still extremize 
St with respect to arbitrary variations of X l (t) and 7Tj(t) with fixed end points 
X|). Doing so one obtains an "almost" extremized action since it is stationary 
for all variations but those of the constant lapse N. 

At this point the "proximity" of 0(i) , p</,(i) with some classical trajectory (f) B (t) , p? (t) 
must be called upon. Indeed, to first order in the change in the lapse around its 
background value N B , one can replace N in the "almost" extremized action by the 
background value N B and obtain 

S T = J q dtfaF+pti- [H G (X\t),n(t)) + H m {(t>,p^X\t))\} (57) 

where X l (t) is evaluated along the stationary orbit from {X\,X^) in a time N B 
which is driven by the off-shell "pressures" dxiH m ((j),p^, X 1 ). The reason why this 
replacement can be performed is the following. Upon searching N, the solution 
of 8nSt = where St(N) is the "almost" extremized action defined above, one 
finds that N is complex function since no classical (i.e. real) solution can be 
found. However, when 0(t) , p^it) is close to (p B (t) , p B (t) , i.e. when if m (0,p^) — 
H m ((fi B ,p B ) < H m ((p B ,p B ) for all t, the real part of N is close to N B and that 
its imaginary part is small. Remember that the location of the saddle is mainly 
determined by the heavy coordinates X % . Thus, one can correctly replace the 
complex by the background value N B . Then, the off-shellness of the matter 
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configuration is made manifest by the fact the H G + H m in eq. fl5?D does not 
vanishes. Not surprisingly, it is through this term that the usual matter action will 
be recovered. 

Having performed this first step, one can indeed proceed to the second and 
third step of Appendix A. The simplest way to proceed consists in "adding" 
H G (X lB , irf) + H m ((f) B , p? , X B ) = to the integrand of eq. fl57D and implementing 
directly eqs. fl5Tj ) and fl52"| ) in the resulting action. (Both equations rely only on 
the stationary character of the orbits X l (t) and X lB (t).) After some elementary 



algebra, one recovers eq. fl54j) 



What we have learned from this generalization is that the determination of the 
linear gravitational response does not require on mass-shell matter sources. Thus 
the matter action which determines the phase of off-shell configurations in a path 
integral formulation of quantum matter propagation in a given background can be 
considered as the action emerging, through a first order variation, from a more 
quantum framework in which gravity responded to these off-shell configurations. 
Even though the imaginary part of N is negligible for heavy coordinates X % , strictly 
speaking this imaginary part is present. So what does it mean? Simply that in 
the path integral of quantum gravity, these gravito-matter configurations will be 
exponentially reduced by this imaginary contribution with respect to the BFA 
formalism in which one has discarded this contribution. 
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